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Abstract 



x 

We study the behaviour of scalar fields on background geometries which undergo quantum 
tunneling. The two examples considered are a moving mirror in flat space which tunnels 
through a potential barrier, and a false vacuum bubble which tunnels to form a black 
hole. WKB approximations to the Schrodinger and Wheeler-DeWitt equations are made, 
leading one to solve field equations on the Euclidean metric solution interpolating between 
the classically allowed geometries. The state of the field after tunneling can then be 
determined using the method of non-unitary Bogolubov transformations developed by 
Rubakov. It is shown that the effect of the tunneling is to damp any excitations initially 
present, and, in the case of the black hole, that the behaviour of fields on the Euclidean 
Kruskal manifold ensures that the late time radiation will be thermal at the Hawking 
temperature. 



1. Introduction 



The radiation of particles from matter evolving along a classical trajectory has been 
heavily studied in recent years. Less well studied is the radiation accompanying quantum 
tunneling from one classically allowed trajectory to another. The following question is 
of interest: if a matter system impinges upon a potential barrier with a radiation field 
in a certain state, what is the state of the field given that the matter is subsequently 
observed to be on the other side of the barrier? A method to answer this question in 
the context of false vacuum decay in flat space was developed by Rubakov Q and has 
been generalized to include gravity as well as topology changing processes |2], |3|. The 
spectrum of radiation is found by solving an imaginary time Schrodinger equation, the 
occurrence of which leads to novel features. Instead of solving field equations in real 
time, one is naturally led to consider propagation on the Euclidean solution interpolating 
between the two classical trajectories. As phase factors in real time are converted into 
exponential damping factors in imaginary time, the resulting particle creation can be 
distinctly different and is accompanied by the systematic supression of excitations present 
before tunneling. 

Given this situation, it is natural to ask how the radiation from black holes might 
be affected by the presence of tunneling. If we consider a distribution of matter, initally 
outside its Schwarzschild radius, which tunnels through a potential barrier to form a black 
hole, the conventional calculation ||] of the radiation does not apply. On the other hand, 
it would be shocking if the same answer was not obtained for the radiation at late times, 
as this is thought to depend only on the hole's late time geometry and not on its history 
at early times. Here we compute the radiation for this process and show that while the 
Euclidean time evolution has an effect at early times, it has none at late times so that 
the standard result is in fact obtained. 

In order illustrate the technique of Ref. || in a simpler setting, we first study the 
effect of tunneling on another well known radiating system — the moving mirror [^]. We 
show in Sect. (2) how an imaginary time Schrodinger equation emerges from a Born- 
Oppenheimer approximation, and use this result to calculate the shift in the spectrum 
of radiated particles as a result of the tunneling. It is shown that the initial spectrum is 
shifted to favor low energy excitations, as is understood by realizing that the probability 
to tunnel is increased if energy is transferred from the radiation to the mirror. 

In Sect. (3) this approach is extended to include gravity in asymptotically flat space. 
A WKB approximation to the Wheeler-DeWitt equation, as considered in ^ 0|, is used 
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to obtain an imaginary time Schrodinger equation which can then be solved as before. In 
Sect. (4) we use this result to examine the radiation from a black hole which is formed 
by tunneling. In particular, we consider the tunneling of a false vacuum bubble, a system 
extensively studied in Refs. [|J This example involves a complication due to the 

peculiar structure that arises; Refs. [13, 14] show that the sequence of three-geometries 
encountered during tunneling can not be stacked together to form a manifold. Employing 
a slight modification of the standard approach, we show how the behaviour of fields on 
the Euclidean Schwarzschild manifold protects the late time radiation from being affected 
by tunneling. An intuitive reason for this is that the bubble's tunneling probability is 
unchanged by the presence of Hawking radiation, which involves the creation of pairs of 
particles with zero total energy. 



2. Tunneling Mirror 



Consider a mirror moving in a one dimensional potential in the presence of a massless 
scalar field. The Schrodinger equation for this system is 
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(2.2) 
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5(f)(x) 2 

Note that \1/ is a function of the mirror coordinate x m , and a functional of the field 
configuration <p{x). The mirror boundary condition is imposed by demanding that the 
field vanish at x m , 

*[(f>,x m ;t} = if (f>{x m ) ^0. (2.4) 

The system is solved by assuming that the backreaction of the field on the mirror is a 
small perturbation of the mirror's motion, and that the mass and momenta of the mirror 
are large enough that it can be described by a well localized wave packet. In this domain 
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the system admits a Born-Oppenheimer approximation, which amounts to an expansion 
in 1/m. In particular, we seek a solution to the time independent Schrddinger equation 

[H m + HjMt, x m ] = E^![(j>, x m ] (2.5) 

valid to zeroth order in 1/m. Following Refs. [|T], |5| the Born-Oppenheimer approximation 
is implemented by writing \l/ in the form 

x m ] = ^vv(x m ) e lSM x m ] (2.6) 

where ipvv is a slowly varying function to be identified with the Van Vleck determinant. 
To lowest order in 1/m, (|2.5|) reduces to the Hamilton- Jacobi equation. 



dS V 



2m ydx m 

since dS/dx m , V(x m ) and E are all of order m. 
To zeroth order: 



+ V(x m ) = E (2.7) 



i d 2 S i dS dipvv r , , , s x 

77- XT ^ X 0> "J j XVP, X m\ (2.8) 
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i dS d 
ifov-j — « — %[</>, arm] + Vvv a™] = 0. 



i/jvv is chosen so that the first two terms cancel, leaving 



i dS d 

Htp x[<f>, %m] = —i — ~ — x[0, x m ]. (2.9) 



This can be put in a familiar form by defining the time variable r(x m ). In a classically 
allowed region, where E — V(x m ) > and dS/dx m is real, r is defined by 

dr m 

allowed regions (2.10) 



doCfn dS I doc 

whereas in a classically forbidden region with dS/dx m imaginary. 

drE . m 



ddbyyi dS j doc 
The resulting zeroth order equations for are: 



i— — — — forbidden regions. (2-11) 



d 

x[4>, t] = i-g- x[</>7 T ] allowed regions (2-12) 



_d_ 

dTp 



x[<fi, t~e] forbidden regions. 



(2.13) 



These are the fundamental equations governing the evolution of the scalar field in the 
presence of the mirror. In the allowed regions we have recovered the time-dependent 
Schrodinger equation with the postion of the mirror playing the role of a clock, whereas 
in the forbidden regions we have obtained a diffusion equation, which we shall refer to as 
the Euclidean Schrodinger equation, with the Euclidean time te measuring the position 
of the mirror in the potential barrier. 

Now, choose the potential to be of the form illustrated in Fig. 1 and let the mirror 
come from right to left. In the allowed region to the right of x l m the state r] obeys 
the normal Schrodinger equation, and so standard methods can be used to find 
Between x l m and x^ m the mirror is in a forbidden region, so the state evolves according to 



1 



dx 
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x[<t>, t e] 



(2.14) 



with x[0) r j?] = x[0) rl ]- We wish to solve this equation in order to find the state at 
the final turning point x' m . It is useful to transform the mirror to rest by defining the 
coordinate 

y(x, t e ) = x -x m (r E ) (2.15) 
in terms of which the Euclidean Schrodinger equation is 
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X[<f>, M 



X [<P,te] (2.16) 



;2.17) 



The solution is 



x[<P,Te] = Texp 



H${T E )dT E 



x[<P, 



Ti 



t^(TE,4)x[^,4] 



:2.18) 



Here T represents time ordering with respect to r' E . The crucial point is that the Eu- 
clidean time evolution operator, Ue, is non-unitary. This is natural since we know that 
wavefunctions decay exponentially during tunneling. If Ue was unitary, the easiest way 
to calculate it would be to transform to the Heisenberg picture, solve the field equa- 
tions mode by mode, and compute Bogolubov coefficients. However, as emphasized in 
Ref. M the non-unitarity of Ue implies that the Schrodinger and Heisenberg pictures are 
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inequivalent, making the standard method inapplicable. Instead, one can use the method 
developed in Ref. ]I[ which closely resembles the standard one but is more general. We 
first describe the state right before tunneling. For convenience, set x l m = r % = r E = 0. 
Let £ u (x,t) be a complete set of positive norm solutions to the Klein-Gordon equation 
which vanish vanish at the mirror: 



2 



+ 



d 2 



dr 2 dx 2 



£u(x,t) = 



i dx 



Cj(x m (r),T) = 0. 



(2.19) 

(2.20) 
(2.21) 



The set of allowed frequencies uj is taken to be discrete, and J2u represents summation 
over this set. The field operators can then be expanded in terms of these modes: 



x 



with 
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[X,T) 



E 



d d 
au>Q-Cj(x, t) + a>l—£u(x, t) 



(2.22) 
(2.23) 



Now define Euclidean fields <fi E (y,TE), n E (y,TE) which agree with 0(x, r), ^(y, r) at 
r = t e = 0, but evolve according to 



j> E (y, r E ) = Ue\t e , 0) 4> E (y, 0) U E (r E , 0) 
7rf {y, t e ) = U e \t e , 0) rrf {y, 0) U E (r E , 0). 



(2.24) 
(2.25) 



We will calculate U e (t e , 0) by first finding 4> E (y, t e ), Ti E {y,T E ). The field equations for 
these operators are 

^ E ,H E 



3t r 



. a Ej , dx m 



So 
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(2.28) 
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and 
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Equation ( |2.29| ) can be obtained by varying the action 

S = 2 / dydrE ^9e 9e^A^ 

with the Euclidean metric 



ds E = g E dx^dx v = dr E + 2 — — dx dr E + dx 2 . 

dr E 



[2.29) 



[2.30) 



(2.31) 



i} E can be expanded in terms of modes f w which satisfy the Euclidean Klein-Gordon 



equation ( 2.29|) and which vanish at y = 0, 

4> E (y, t e ) = t e) 



(2.32) 



Or, 



d dx d 

fu(y, te) - -j—^-foj(y, t e) 



dr E dy' 



;2.33) 



As the Euclidean Klein-Gordon equation is elliptic, one cannot in general impose 
Cauchy boundary conditions at t e = on f u . The resulting solutions would not satisfy the 
mirror boundary condition. With the appropriate boundary conditions, either Dirichlet 
or Neumann, imposed at t e = and t e = t e , a detailed calculation is, of course, required 
to find fu for a generic mirror trajectory. We shall take the solutions as given and only 
use their specific forms in a region far from the mirror, where they are simple. 

Now, using the condition that the two sets of operators 0, ix^ and <fi E , n E are equal 
at r = t e = 0, and taking inner products, the operators b w can be expressed as a linear 



combination of d U) a* : 



(2.34) 



Then using 



4> E (y,4) = U E \r E ,0)My,0)U E (T E ,0) = U E \r f E , 0) 4>(y, 0) U E (r E , 0) (2.35) 



and the analogous expression for j} E , the following equations for U E are obtained: 
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and 



E 



£ [t^rf, 0) a* ^(r^, o) Uv, o) + E^Te, o) at U E (r f E , o) o) 



d 

ue'H, o) a* u E (4, o) — ^(y, o) + u- E \rl o) a+ ^(r|, o) —C{y, o) 



Again taking inner products, this leads to relations of the form 



Ue X {te, 0) a* C/e(t|, 0) = ^ [itwaw' + 
Then it can be shown that [DJ1 



U e (t e ,0) = const. x:exp^J^ 



-D^ajAj + F^ia^hrji + -G^a^a 



1 



u;u; ^o; a; 



' war u, w u, W 



where the matrices .D, F, and G are defined by 

D = vz- 1 ; F=(z T )~ 1 -l ■ G = -z- 1 w. 
The state after tunneling is then determined, 

X (4)) = u E (4) |x(o)} 



(2.36) 



(2.37) 

(2.38) 
(2.39) 

(2.40) 

(2.41) 
(2.42) 



and is expressed in terms of occupation numbers with respect to the modes ^(y, 0), where 
now y = x — x* m . All of the information about the final state is contained in the matrices 
D, F, and G, which are in turn given in terms of inner products between the modes f w 
and £ w . 

As a simple application of these formulae we will calculate the shift in the spectrum of 
outgoing particles which are far from the mirror at the time of tunneling. It is assumed 
that the mirror was initially at rest and the field in its ground state. The mirror subse- 
quently accelerates in the potential V(x m ) until it reaches the classical turning point x l m . 
It is well known that as a result of the mirror's acceleration, a flux of outgoing particles 
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is created whose spectrum is calculable by standard methods f7j. Outgoing particles far 
from the mirror are wavepackets composed of superpositions of plane waves, 



1 



2v^ 



-iu>(r~-x) 



2.43) 



The spectrum of outgoing particles located at x = x ^> uj 1 atr = 0is written as 



E^(K»iK}) 



(2.44) 



where {n^} is a set of occupation numbers and S s ({n u }) is the amplitude for the set to 
occur. 

Far from the mirror, the modes f u are easy to calculate since the mirror boundary 
condition is irrelevant. They are of two types, 



f 



+ 



o L0T E -\-iiJX 



1 



3 wte +iu) (y+x m (t e ) ) 



Then (ft, ft and <p , are equal at r = r E = if 



(2.45) 



K = ow ; &j~ = oj,- 



;2.46) 



Equation (|2.38|) gives 



leading to 



and 



U E \r E , 0) a^U E (4), 0) = e-< +i ^' 
U E 1 (r E} 0)aiU E (T E ,0) = e^ 



D = G = ; F w 



Ue(t e ,0) = const, x : exp^ 



const, x : exp y j 



-iux f m 



a^a^ :: exp^ 



e _WT « - 1 



;2.47) 



'2.48) 



;2.49) 



The first factor is a translation operator which expresses the state in terms of the modes 
^(x,0) instead of £ w (a; + x^,0), and the second factor acts on a state to give 
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e E ({ n <»}) T E \{n u }), where E ({n^}) = J2 n u& is the energy of the state. Therefore, the 
state after tunneling is 

const, x ]T e~ E ^ T e S s (K» \{n„}) . (2.50) 

{riu,} 

The result of the tunneling is simply to shift the spectrum from S x to e~ E( - n ^ T E S x . 

It is not difficult to understand this result. Since the total energy is fixed , the 
state before tunneling is given by a superpostion of the various ways of distributing 
the energy between the mirror and the radiation. As the mirror's probability to tunnel 
depends exponentially on its energy, we expect an inverse exponential correlation between 
tunneling and energy in radiation. Thus an observer measuring the spectrum of radiation, 
conditional on the mirror tunneling, finds the result ( |2.50| ). Far from the mirror the shift in 
the spectrum depends only on r^, the amount of Euclidean time spent during tunneling. 
This is because the tunneling amplitude in the WKB approximation is e~ s , and the 
derivative of S with respect to energy is just the Euclidean time. 

If we were to identify the Euclidean time with an inverse temperature, the shift would 
become a Boltzmann factor. This makes it easy to generate thermal distributions of 
radiation. Specifically, if the distribution before tunneling was a constant, then after 
tunneling tracing over the states of the mirror would yield a thermal density matrix for 
the radiation. A number of authors have been led by this fact to seek a connection 
between the thermal radiation that arises in the contexts of cosmology and black holes 



and an occurrence of tunneling |3], [18], Such a connection relies upon assumptions 
about what is on the other side of the barrier and what the spectrum of radiation is there. 
In this work we only consider situations where there is a well defined classical trajectory 
on either side of the barrier; we are interested in the case in which there is collapsing 
matter on side of the barrier and a black hole on the other. The treatment of this process 
requires an extension of the previous method to include gravity. 



3. Application to Gravity 

In this section we make a WKB approximation to gravity in a manner which directly 
parallels that for the moving mirror. The starting point for the canonical quantization of 
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gravity is to write the metric as 

ds 2 = - (N^tf + h i:j (dx l + N l dt\ (dx j + N j dt) . (3.1) 
With this definition, the action for gravity plus matter takes the form 

2 

S = — — / d A x^f^g (1Z — 2 A) + Sm + boundary terms 

167T J 

= J d 4 x (ir^ft + ^jh ij - N^t - NiT-lf) + boundary terms (3.2) 



with 



Ht = ^/T* {h ik h ol + huh jk - hijhu) <K ij <K M - ^h* ( 3 TZ - 2 A) + H tM 

' (m G ijM ^ kl - ^ ("Jl - 2A) + H tM (3.3) 



m~ 1d7T 

n = -27r%. + H M (3.4) 

Covariant differentiation and the raising and lowering of indices are performed with respect 
to the spatial metric h^. <pi refer to arbitrary matter fields. 

The boundary term in the action is determined by requiring that no such term arise 
in the variation of the Hamiltonian [15|. Restricting to asymptotically flat metrics with 
N f — > 1, Ni — > as r — > oo, but allowing for time translations at spatial infinity, the 
appropriate boundary term is numerically equal to the mass of the system. Thus the 
Hamiltonian is 

H = Jd 3 x [N*7U + Nirf] + M. (3.5) 
Varying the action with respect to AT* and Aj yields the constraints 

H t =H i = 0. (3.6) 

The system is quantized by making the replacements 



7T 



~ l TT~ 5 n <h ~^ ~ l TT (3.7) 



' Shij 



and demanding that the state satisfy the constraints Tit^ = Ti, 1 ^ = as well as the 
Schrodinger equation H^> = id^f/dt. 
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The constraints 



2i 



Sh 



n 



M 



ij/ \j 



* = 



(3.8) 



enforce invariance of the state under spatial reparameterizations, and the constraint 



167r 5 5 
m 2 13 kl 5h {j 5h kl 



m p - 



167T 



hh ( 3 n - 2A) + n tM 



* = o 



(3.9) 



is the Wheeler-DeWitt equation. Proceeding as before, we seek a semiclassical solution 
of the form 

* [hij, <j>i\ = ipwihij] e im ? s ^ X • 
At first order the Einstein-Hamilton- Jacobi equation is obtained: 

,2 



167r^ 5S 5S m v x ,„ „ , 



2 J «3*i^.. ^ 



167T 



0. 



Zeroth order yields 



provided ?/Vv satisfies 



167T 5x . V/ n 



771 



V kl XU XU WW ^ijkl~ — U. 



' 6 6 h i;) 



The momentum constraints at first order are 

( SS \ 



'8 h^ Sh H 



5h 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



and at zeroth order are 



2i 



6h 



+ n i M x = o- 



13 / \j 



(3.15) 



Equations (p. 12 ) and ( 3.15 ) describe how the matter wave function evolves as the 
spatial geometry changes. Quantum field theory in curved space can be recovered by 
writing x's dependence on h^ in terms of a time functional t[x; hij], and by reintroducing 
a lapse N T and shift Ni, demanding that they obey 



G 



5S 



ijkl 



5hij 16nN 7 



dy 



5h 



5r[y; h ab ] 



(3.16) 
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Then 



,16tt 



N T G ijkl ^--^ + 2iN i 
oha oh 



ij <J'l<kl 



Shin 



—I 



St Shu 



so that the equation for \ becomes 



d 3 x 



tM 



NjTt 



d 



(3.17) 



(3.18) 



The condition ( |3.16| ) agrees with the classical relation between and hij, demonstrating 
that t[x; hy] is the classical time and that ( |3.18[) is the Schrodinger picture version of 
quantum field theory in curved space. 

As with the mirror example, r becomes imaginary during tunneling so we define a 
Euclidean time te along with a Euclidean lapse N TE = iN T , in terms of which x obeys 



- / dx[N^H tM + iNi# M ] X [<Pi, r E ] 



_d_ 

dTp 



-XW,te\ 



For a massless scalar field with action 



we have 



Tit, 



i (h-l 

rii 



(3.19) 

(3.20) 

(3.21) 
(3.22) 



To evolve \ through the tunneling region one is required to calculate the Euclidean time 
evolution operator 



Ue(4,t e ) 



with 



m 



Texp 



x S 2 



j 

E -~,E 



KdT E 



d 3 x [N TE (-^ — + ^hWdrfdrf) + Nidvr . 

2 d(p z 2 o<p 



(3.23) 



(3.24) 



As before, one proceeds by defining Euclidean fields obeying ( ^24 , ^25) . In the present 
case the resulting field equations are: 



9e 9E d n 







?rf = i 



N t e \dr E 



- N%c 



(3.25) 
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with 

dsl = gf v dx^dx v = (N TE dr E ) 2 + (dx l + N^te) (dx j + N j dr E ) . (3.26) 

The evolution operator, and therefore the state after tunneling, is determined by solving 
the field equations mode by mode, and repeating the steps leading from ( |2.32| ) to fl2.42| ). 



4. Black Hole Radiation in the Presence of Tunneling 

We can now apply this method to determine how the radiation from a black hole is 
affected by tunneling. It is well known that a black hole formed classically from collapsing 
matter radiates in a complicated manner at early times due to the time dependent geome- 
try, but at late times will inevitably radiate as a black body at the Hawking temperature. 
Is this scenario altered if the black hole is formed while tunneling? We shall show that 
it is not. The form of the late time radiation is insensitive to the hole's unconventional 
history in a way that is consistent with the intuitive picture of Hawking radiation being 
caused by pair production near the horizon. 

We consider the behaviour of a scalar field on the background of a false vacuum bubble 
which tunnels leading to the formation of a black hole. The action for a false vacuum 
bubble in the thin wall approximation is 

S = ^~ fd'x^Fn-^- [ d *xV=9--^[ d 3 A (4.1) 

lu7T J QTT J bubble AlT J wall 

where A/ is the cosmological constant of the false vacuum, and fi is the energy density 
of the bubble wall. The classical solutions for this action have been derived in Refs. IHI- 



121 . In what follows we refer to the treatment of Ref. ||12|| . The spherically symmetric 
solutions are characterized by three parameters: A/, /z, and the total mass M. In addition, 
for given Aj and /i there is a critical mass M cr below which there are two solutions: type 
(a), where the bubble emerges from a singularity with zero radius, subsequently expands 
to a maximum radius, and then recollapses; type (b), where the bubble initially collapses 
from infinite radius, reaches a minimum radius, and then reexpands. Using the results of 
Refs. ]I3|, [L4| , we focus on an expanding solution of type (a) which tunnels to an expanding 



solution of type (b). We confine our interest to the region outside the bubble where the 
metric, written in terms of Schwarzschild time t and r* = r + 2Mln(r/2M — 1), is 

2M\ 



ds< 



(l _ f^j ^dt 2 + drt) + r 2 dn 2 . (4.2) 
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As t and r* cover only part of the complete manifold, we introduce Kruskal-Szekeres 
coordinates, 

ds 2 = {-dT 2 + dX 2 ) + r 2 dn 2 . (4.3) 

r 

The two sets of coordinates are related by 

e r/2M =X 2_ T 2 



\2M 

_ j 4Mtanh- 1 (T/X) if \T/X\ < 1 

\ 4Mtanh- 1 (X/T) if \T/X\ > 1 1 ' 

Using these cordinates the type (a) and (b) solutions of interest are depicted in Fig. 2. 

The tunneling amplitude for this process has been computed by two different methods. 
In Ref. |14[] the solution to the Wheeler-DeWitt equation is found in the WKB approxima- 



tion by solving the Einstein-Hamilton- Jacobi equation ( 3.11 ). Since the solution behaves 
as e~ s , and the tunneling amplitude is given by the ratio of the wavefunction evaluated 
at the initial and final geometries, the tunneling amplitude is 

exp (S[/4f ial ] - S[hf™ 1 ]) (4.5) 

No difficulties arise in this approach; the calculation of tunneling amplitude proceeds in 
a straightforward fashion. 

In Ref. the calculation is performed using the functional integral. In this formalism 
one looks for a manifold which interpolates between the initial and final surfaces and which 
is a solution to the Euclidean Einstein equations. The tunneling amplitude is e~ s , where 
S is the action of the solution. It is found, however, that solving the field equations leads 
to a sequence of three geometries which do not form a manifold. To see this, first note that 
the geometry outside the bubble is Euclidean Schwarzschild space, obtained by t —>■ tig, 
T -> iT E , 

ds l = (i - ^ (dt 2 E + dr 2 ) + r 2 dQ 2 = 32M3 e~ r/2M ( rfT 2 + dX 2^ + r 2 dQ 2 (4 6) 
with 

{Wl ~ ^ = x2 + T e^e = 4Mtan- 1 (7VX). (4.7) 

It remains to describe the motion of the bubble wall. Solving the equations of motion 
leads to the trajectory in Fig. 3. It is seen that the bubble wall crosses the initial surface 
during the course of its motion, creating a situation in which it is impossible to identify 
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a region which is swept out by the evolving hypersurface. Some regions of the manifold 
are crossed twice by the hypersurface, some once, and some not at all. The authors of 



Ref. JL3J call this object a pseudomanifold and give a prescription to calculate its action 
by assigning covering numbers to the various regions, but this is not needed for what 
follows. 

With these results in hand, the technique of Sect. (3) can be used to calculate the 
state of the scalar field after tunneling. It was seen that once the solution of the Einstein- 
Hamilton- Jacobi equation is given, the field wave functional x is fully determined by 
equations ( |3.12p and Q3.15| ). Since S[hij] is calculated in Ref. |T5 |, we have all that we need 



to find x- This would, however, require finding the solution to an unfamiliar functional 
differential equation. To cast it in in the form of the Schrodinger equation a lapse N TE , 
shift iVj and time te were reintroduced leading to the appearance of the Euclidean metric 
gz v . In the present case there is no true interpolating Euclidean manifold, so that any 
choice of N TE and iVj which define a well behaved g^ v will lead to a bubble trajectory 
that is a multivalued function of time. Alternatively, a choice of time functional which 
gives a single valued bubble trajectory will necessarily lead to a Euclidean metric with 
vanishing determinant at some point. In either case, it is not clear that the resulting 
Schrodinger equation is well defined. This is apparent from Fig. 3, where it can be seen 
that boundary conditions imposed on the initial surface and on the bubble wall may 
contradict each other. These difficulties arise as a result of trying to compute the final 
state of the field in one step, which requires a Euclidean manifold interpolating all the 
way from the initial surface to the final surface, and can be avoided by calculating the 
state on a series of intermediate hypersurfaces. In this approach, it does not matter that 
the bubble wall eventually crosses the initial surface since once the state is calculated at 
some intermediate point we can forget about what preceded it. 

For simplicity, we will consider only the s-wave component of the scalar field and 
frequencies high enough such that the geometrical optics approximation is valid. This 
means that the field equation is taken to be 

ggd^t = 0. (4.8) 

The state of the field on the initial surface, t = T = 0, is most conveniently expressed in 
terms of the coordinates r* and t. We divide the modes into ingoing and outgoing, 

e(r*,t) = cu- Mi+r * } 

C n \n,t)=C UJ e-^ t -^ (4.9) 
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and write the field operator as 

#(r„ f)=£ + *) + in - out] • ( 4 - 10 ) 

C w are normalization constants whose values will not be important. We shall only consider 
the in modes as the treatment of the out modes is exactly the same. We also suppress 
the in superscript. 

In the first stage of the evolution the hypersurface is pivoted around r* = by 180°, 
where r b is the position of the bubble wall on the initial surface. The solutions to the 
Euclidean field equations are most conveniently obtained by choosing Cauchy boundary 
conditions on the initial surface, (clearly a valid procedure in this case) 

ti(r*,0) = Ur*,0) ; ^L + (r«,0) = -*-e(r,,0) 

/ w -(r.,0) = C(r.,0) ; ^/J(r*,0) = -i|e(r*,0) (4.11) 

It is also easiest to use the X, T coordinates as they are well behaved everywhere. Since 
the evolution of the hypersurface is simply a reflection about the point X = X b , a mode 
which has the form f(X,T E ) on the initial surface has the form f(—X + 2X b ,T) on the 
new surface. Using the relations 



r, = 4Mln^X 2 + T| ; t E = AM tajT 1 (T E / X) (4.12) 

and that 

f±(r.,t E ) = C u e tb *s- iur - (4.13) 
near the initial surface, one sees that near the new surface, 

ft(X,T E ) = aexp ( j^^ - UMu\n{-X + 2X 6 )) . (4.14) 

Since on the new surface, /+ = (/")* and df+/dt E = —(df~/dt E )*, the evolution oper- 
ator U E is unitary. This means that the state on the new surface has the same form as it 
did on the initial surface, but is now expressed in terms of the modes 

UX, T) = C„ exp ( Zx + 2X* + 4tMUJ H ~ X + 2Xh) ) ■ (415) 
These modes can be approximated near T = as 

. _ J C^~^ if |X|»X b 
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Now it is useful to express the state in terms of modes which are nonzero only inside or 
outside the horizon, 

' D^e^-^ if X < 
if X > 



„>-} " if x < 




A fundamental result ||, [16| in the derivation of black hole radiance is that the vacuum 



state with respect to modes which have a time dependence e tujT is the state 

const, x ]T e -H{n»})/2T H ){rU)< | {?U)> (41g) 

{n u } 

with respect to the modes rj^ and rj^. The sum runs over all sets of occupation numbers, 
E = X)n w tJ, and Tjj = l/8nM is the Hawking temperature. Further, near the horizon, 
any deviation of \x) from the vacuum state can be ignored because of the arbitrarily large 
redshift as r„ — > — oo. Far from the horizon £ w and rj< agree so the form of the state is 
unchanged there. 

Now the hypersurface can be evolved the remainder of the way. If we restrict our 
attention to the region X < X b , then the motion of the hypersurface is simply a trans- 
lation, tE — ► tE — AtE- This causes states with time dependence e lwt to be damped by a 
factor e~ wAii5 , and states with time dependence e~ luJt to be amplified by a factor e wA * B . 
Near the horizon, the state \x) consists of pairs of positive and negative frequency states 
according to (|4.18 ). One member of the pair is damped but the other is amplified by a 



compensating amount so as to leave the state |%) unchanged. The final state of the field 
can then be summarized as follows. Far from the hole, where there is no pairing, the 
initial state is damped: 

E S(K)) IK)) — > const - x E e' E ^ At - |K}) . (4.19) 

Near the horizon the final state is given by fl4.18p . This is true for both the in and out 
modes, so an observer stationed on either side of the horizon would observe a thermal 
distribution of both ingoing and outgoing particles. As time passes, all of the ingoing 
particles will eventually cross the horizon and be swallowed by the hole, whereas the 
outgoing particles will propagate out to infinity where they can be detected at arbitarily 
late times as a flux of thermal radiation at the Hawking temperature. 
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5. Conclusion 



It was shown that the standard picture of black hole radiance is unchanged by tun- 
neling. At late times, the hole radiates just as it would have had it been formed from a 
classical collapse. This makes sense if one thinks of Hawking radiation as pair production. 
The probability of tunneling is not affected by the creation of a pair, since the pair has 
zero total energy. From this point of view it is also clear that what happens at early times 
cannot possibly affect the late time radiation, since the produced pairs only see the late 
time geometry. The conventional derivation of radiance obscures this point somewhat and 
it seems desirable to find an approach which makes this feature manifest from the outset. 
For the two systems considered here, and presumably this is true in general, the effect 
of the tunneling was to shift the distribution of any particles that were present before 
tunneling. In the present case initial excitations were damped because the final surface is 
rotated clockwise relative to the initial surface. A counterclockwise rotation would have 
led to amplification. In [13 numerical investigations are quoted which show that the 
rotation is always clockwise for the false vacuum bubble. One is led to speculate whether 
this is a general phenomenon — whether all tunneling transitions lead to damping. 
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Figure Captions 

1. A generic mirror potential. The turning points for energy E are indicated 

2. The type (a) and (b) solutions. The heavy lines represent the bubble trajectory, and 
the dashed lines are the initial and final surfaces of the tunneling solution. In these figures, 
only the regions to the right of the trajectory are of interest, as they are outside of the 
bubble. 

3. Bubble trajectory in Euclidean Schwarzschild space. 
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